Precision cosmology with a combination of 
wide and deep Sunyaev-Zeldovich cluster surveys 
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We show the advantages of a wedding cake design for Sunyaev-Zel'dovich cluster surveys. We 
show that by dividing up a cluster survey into a wide and a deep survey, one can essentially recover 
the cosmological information that would be diluted in a single survey of the same duration due 
to the uncertainties in our understanding of cluster physics. The parameter degeneracy directions 
of the deep and wide surveys are slightly different, and combining them breaks these degeneracies 
effectively. A variable depth survey with a few thousand clusters is as effective at constraining 
cosmological parameters as a single depth survey with a much larger cluster sample. 

PACS numbers: 



Introduction. — The abundance and redshift distribu- 
tion ^ of clusters are potentially important probes of 
cosmological parameters [lH5|. Recently, number counts 
from cluster surveys in optical (e.g RCS [3] and SDSS 3) 
and microwave [9| have been used to constrain cosmol- 
ogy. In the coming decade, dedicated surveys in differ- 
ent wavebands (for example, ACT, DES, eROSITA, etc.) 
plan to use clusters as probes of precision cosmology. A 
correct cosmological interpretation of the observed ^ 
depends on precise knowledge of the limiting mass of the 
survey at a given redshift. Unfortunately, the mass of 
a cluster is not directly observable, but has to be ob- 
tained through proxy observables such as X-ray surface 
brightness and temperature [lO|, Sunyaev Zeldovich (SZ) 
decrement ll|, ll2| , cluster richness 13|, |lj] and lensing 
ISl Il6| . Uncertainties in the observable to mass conver- 
sion create degeneracies between cosmology and clus- 
ter physics and degrade cosmological constraints. There 
have been a number of attempts in the literature to break 
such degeneracies, for exanxple, through the so called 
'self-calibration' techniques [a, l6lll7|. Other approaches 
include an 'unbiased' mass follow-up of a subsample of 
the survey clusters IT], |l8| or better theoretical modeling 
of clusters to predict the form of mass-observable scaling 
relation 19I421J. All of these are equivalent to putting 
priors on the parameter space of the mass-observable re- 
lations. One can also try to optimize the cluster surveys 
so as to get the best possible survey yield |22l|. 

In this Letter, we show that by simply dividing the 
total observation time between deep and wide surveys, 
parameter degeneracies can be effectively broken leading 
to tight cosmological constraints. This is without the 
need of any external information or a costly mass follow- 
up program. Throughout this paper, we assume a spa- 
tially flat cosmological model with a constant equation 
of state w for dark energy with the fiducial cosmological 
parameters taken as h = 0.717, fim — 0.258, w = —1, 
fib — 0.044, Us = 0.963 and as — 0.796 consistent with 



the WMAP 5-year results. 

Preliminaries. — The redshift distribution of de- 
tectable clusters in a survey over an area Afi is. 
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is the comoving volume element and -^ is 
the cluster mass function (taken from simulations by 
Jenkins et al. [23[). A mass limit Miim(z) is effec- 
tively imposed through the complementary error func- 
tion f{M\Mii,j^{z)) [2J|, which is used to model a 20% 
logarithmic scatter in the mass-proxy relation centered 
about Miim. Figure [1] shows the plots of Mum ver- 
sus z for various array sensitivities of the ACT survey 
along with the corresponding redshift distribution of clus- 
ters. For an SZ survey, the flux-mass scaling relation 
is used to determine Mii„^{z) fig, fluxsz(z,r^)d^(z) = 
2.699 X 10^^g{v)AszM^aoE^/^{z) (1 + zy , where g{iy) is 
the frequency dependence of the SZ distortion, Af2oo is 
in units of Mq and dA is the angular diameter distance 
in units of Mpc. The Hubble expansion is parameterized 
as H{z) = HqE[z). The parameter 7 captures deviation 
in the redshift evolution of the cluster scaling relation 
from the self similar case, and we assume 7 '^ 0. We 
choose the SZ scaling relations to be compatible with re- 
cent observations 2l|,|25| with log(Asz) = -28.067 and a 
= 1.612. 

Majumdar and Mohr IT], ll8| have shown that a sub- 
stantial improvement in the cosmological parameter con- 
straints can be obtained by adding a follow-up of a 
fraction of the cluster sample to calibrate the mass- 
observable relation through detailed observations in X- 
Ray, SZ or galaxy spectroscopy. We construct a mock 
follow-up survey with ^ 100 clusters having 'unbiased' 
mass measurements in the range 2 — 10 x 10^"*^-"^ Mq 
between 0.3 < z < 1.1 with 20-40% error on masses. 
We use both Fisher matrix and MCMC techniques, de- 
tails on which can be found in many references, e.g. 




FIG. 1: The cluster redshift distribution ^ (thick lines) and 
Afiim (thin lines) vs redshift, for a SZ survey with array sensi- 
tivities of 3.54 (solid), 10.61 (dot-dashed), 21.21 (dotted) and 
34.12 (dashed) ^K^s. 



P, i, 0, [II m, [13. Full MCMC runs are made to esti- 
mate the covariance matrix between cosmological param- 
eters and {Asz,ci,"f) from follow-up. The full parame- 
ters space is spanned by {h, Ug, ^h, ^m, w, o'SjAsz, a, 7) 
and we use Gaussian priors with standard deviations of 
0.027, 0.015 and 0.003 on the parameters h, Ug and ilb 
respectively. 

We assume the survey instrument to be a telescope 
with an angular resolution of 0fwhm with a bolomet- 



ric detector array of equivalent array-sensitivity (noise- 
equivalent temperature) aa- For a survey of area A and 
duration iobs, the time spent on each of the Npi^ = 



^ ^/^fwhm sky-pixels is tpi^ = iobs/^pix, giving a noise 



per sky-pixel of At — cra/\/^pix- In particular, we as- 
sume a survey of area A = 2000 sq. deg. at 150 GHz 
with 6'fwhm = 1' and a duration of 10^ seconds. We con- 
sider three cases for the array-sensitivity — 3.54, 10.61 
and 21.21 /iK^/s. In perspective, for its 2008 season, 
the ACT had an array sensitivity of ^ 34 iiKy/s at 148 
GHz. These array sensitivities would correspond to flux 
limits of 2.7, 8.3, 16.6 and 26.7 mJy respectively for a 
single area of 2000 sq. deg. in the ACT survey. The 
SPT has similar configurations. Given a thermal noise 
At per sky-pixel, the flux limit is computed as follows: 
we assume that the SZ decrement/increment has to be 
integrated over roughly 30 pixels to detect a cluster, so 
the 5ct flux limit becomes 5-\/30At converted from fiK 
into mJy through the derivative of the blackbody func- 
tion [28|. However, flux sensitivity may be limited by 
the covariant noise from CMB as well as confusion aris- 
ing from point sources; we use a fixed (irrespective of how 
deep we go on a patch) limiting flux value of 2.5 mJy as 
an absolute lower cutoff. 

Results. — 

We illustrate the strength of a wedding cake survey 
by breaking the fiducial 2000 deg^ survey area into two 
patches keeping the total observing time fixed. We use 
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FIG. 2: Parameter degeneracy in the ag-'y plane for a wedding 
cake survey with /time = 0.75 and r = 1/9. The thick ellipses 
indicate constraints when mass follow-up is added, while the 
thin cigar-like ellipses are the constraints for unknown cluster 
evolution. For these, the dashed lines are the constraints from 
the small area survey, while the solid lines are constraints from 
the large area survey. The light (dark) filled ellipse shows the 
combined constraint from the two surveys with mass follow- 
up (unknown evolution). Note that, the intersection of the 
two cigars give the tighter constraints. 
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FIG. 3: Constraints on w (upper panel), and ag (lower panel) 
from only number count observations for a two component 
survey design as a function /time . The curves are plotted for - 
CTa(^K^) = 3.53, r = 1/9 (dot-dashed); aa = 10.61, r = 1/9 
(solid); aa^ 21.21, r = 1/9 (dashed), aa = 10.61 & r = 1/4 
(dotted), and aa = 31.83 &i r — 1/9 (thin). 



/time to denote the fraction of total time spent on the 
smaller area, and r to denote the ratio of the smaller 
to the larger survey area. Varying /time amounts to dif- 
ferent integration time on each patch and hence differ- 
ent flux limits. Different flux limits lead to changes in 
the shape of ^ (see Figure [T|), so that the two patches 
yield different parameter degeneracy directions. When 
the constraints are combined, the parameter degenera- 
cies are effectively broken, leading to tight constraints on 
cosmological parameters. This is evident in Fig. [2] which 
illustrates how the notorious cosmology-cluster physics 
degeneracy (such as erg ~ 7) can be broken. A multiple 
depth survey is akin to having surveys with multiple lim- 
iting masses or Mn^. This provides binning in both the 
observable mass proxy as well as redshift. The binning 
in the mass is what helps in breaking the degeneracies, by 
extracting out more information from the survey [2J| . 

a) Wedding cake only. — The marginalized constraints 
on w and as from an SZ survey like the ACT for vari- 
ous detector noise levels as a function of /time are shown 
in Fig. [3l It shows that having a single sky area to a 
fixed depth with no prior knowledge of cluster parameters 
gives no constraints on cosmology, in qualitative agree- 
ment with previous work [17|. However, the constraints 
dramatically improve as soon as the survey is broken up 
into wide and deep parts. For our survey of total dura- 
tion 10'^ s and with cluster yield A^ci ^ (3200, 2700) for 
the wide and deep parts (for /time ^ 0.75, r — 1/9 and 
aa = 3.54 /xK-^s) one can get {AQm, Aw, Acrg) = (0.105, 
0.115, 0.116) just from ^. For a similar exercise, a sin- 




FIG. 4: Constraints on w in a wedding cake survey. The 
main figure shows the constraints when {Asz,a) are known. 
The thick (thin) curves correspond to an unknown (known) 
cluster evolution respectively. The inset displays constraints 
when (y4sz,a,7) are free parameters but mass follow-up is 
added. The curves are plotted for aaifjKy's) — 21.21, r — 
1/9 [dashed); aa = 10.61, r = 1/9 (solid); aa = 10.61, r = 
1/4 (dot-dashed), and cr„ = 31.83 & r = 1/9 (thtn). 



gle depth survey of 2000 deg^ with N^ '^ 7500 would get 
(Aan,Aw,A(T8) = (0.772, 1.41, 3.31)! Thus, moving 
to deep+wide survey acts like having an additional mass 
follow-up (or 'self- calibration') without any extra effort. 
Note that the bump at /time = 0.1 for r = 1/9 corre- 
sponds to having the same flux limits for the two surveys 
such that the parameter degeneracies are the same and 
there is no improvement in the constraints. For a higher 
area ratio, r = 1/4, the bump is shifted to the right, as 
expected. 

b) Breaking degeneracies: Mass follow-up vs wedding 
cake survey. — Our lack of knowledge in the evolution 
in the mass-observable relations can lead to cosmolo gica l 
constraints getting weaker by factors of 4 or more Il8| , 
which can be restored through either a mass follow-up 
program or using cluster bias as a proxy for mass (i.e, 
self-calibration). However, the former requires a costly 
follow-up effort and 'unbiased' mass estimators, while 
for the latter, one needs to be careful about scale de- 
pendent bias. Here, we put forward a simple solution for 
this problem. Following 17|, we parameterize the red- 
shift evolution of mass - observable relation with a factor 
(1 -I- z)^ in the scaling relation. To start, we assume that 
Asz and a are known perfectly. Fig. |4] shows the con- 
straints on w when 7 is known/unkown. We find that a 
two component survey improves the constraints on w by 
a factor of 4 for the case when cluster evolution is un- 
known. This is mainly by breaking the degeneracy exist- 
ing between the parameters 7 and cg. As shown in Table 
m the constraints on w improve from 0.4 to 0.095 when 
we have a two component survey. Note, that the results 
from MCMC (triangles and circles in Fig.|4]) confirm the 
trends seen from the Fisher pipeline. The inset in Fig. |3] 
shows the constraints obtained on the parameter w when 
both cluster structure and evolution is unknown but we 
have some information on all the three cluster parame- 
ters from mass follow-up, the currently preferred method 
for upcoming cluster surveys. The tightest constraints on 
w are obtained when the smaller patch is observed be- 
tween 70% and 95% of the total time, depending of the 
array-sensitivity. Adding follow-up observations to clus- 
ter counts already breaks the degeneracy in cosmology- 
cluster physics, but a wedding cake design can still im- 
prove w constraints by a factor of 1.5 or more. Again, 
we confirm the trend seen in Fisher matrix forecasts with 
MCMC simulations. A detailed comparison of Fisher-vs- 
MCMC for cluster cosmology will be given in [29| . 

Discussion. — We have shown that by breaking up a 
cluster survey into a wide and a deep parts, tight cos- 
mological constraints can be obtained without the need 
for costly mass follow-up observations or self-calibration, 
which may be limited by our understanding of scale- 
dependent bias. For different scenarios — having only 
number counts, having a follow-up survey, or known scal- 
ing relations with unknown redshift evolution, we always 
find an improvement in the constraints on the dark en- 



TABLE I: Comparison of la parameter constraints from a 
cluster survey of array-sensitivity 10.61 ^Ky^ and duration 
lO'^ seconds. The boldface numbers with similar brackets are 
for easy comparison of a single component vs two component 
survey. 



Single area (2000 deg 


^) survey 






unknown 


unknown 


with mass 


Parameter 


scaling 


redshift evolution 


follow-up 


AQn, 


0.676 


0.019 


0.032 


Aw 


1.343 


(0.400) 


[0.147] 


Ao-8 


2.849 


0.121 


0.028 


Wedding cake (1800+200 deg^) survey 




unknown 


unknown 


with mass 


Parameter 


scaling 


redshift evolution 


follow-up 


Ailm 


0.105 


0.009 


0.030 


Aw 


[0.115] 


(0.095) 


0.088 


Ao-g 


0.116 


0.019 


0.028 



ergy equation of state w. Our finding are sumniarized 
in Table HI which clearly shows that a two component 
survey is always superior to a single area survey. For 
an SZ survey like the ACT, just number count obser- 
vations with the two component design are capable of 
constraining cosmology very well (AJ7m = 0.105, Aw — 
0.115 and Aug = 0.116), when a single area approach 
would produce cosmologically uninteresting constraints. 
These results can be further improved by having follow- 
up observations of ~ 100 clusters. A wedding cake survey 
is also able to effectively break the degeneracy between 
non-standard evolution 7 and as, and provide factors of 
several improvements in constraints over a single area 
survey (third column of Table |l| . 

With the advent of a large number of surveys in the 
near future, we should be able to use clusters as indepen- 
dent probes of cosmology. Precision cluster cosmology 
is traditionally thought to require either a good under- 
standing of cluster physics, or detailed mass follow-up 
observations of a fraction of the detected clusters. We 
show that having a wedding cake survey strategy is an 
effective way to beat down cosmological parameter de- 
generacies. Our work provides an important new insight 
into developing future SZ survey designs. 

The authors would like to thank David Spergel for dis- 
cussions. SM also thanks Anya Chaudhuri, Christoph 
Pfrommer and Joe Mohr for many discussions on scaling 
relations over the years. 
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